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We derive the scaling dimension of boundary operators in the AdS/CFT correspondence using 
a path integral treatment of the boundary-to-boundary propagators of their dual fields. We then 
apply the same technique to AdS/QCD where scaling dimensions are anomalous. In particular, 
we compute the two-point correlation function of gluon field strength operators, for which it is 
prerequisite to know the flow of the anomalous scaling dimension. The results are in very good 
agreement with quenched lattice QCD data, thus confirming the functional form of the scaling 
dimension. 



In this paper, we solve two closely related problems. 
First, we extend the formula for the scaling dimension 
of boundary operators known from gauge/gravity duali- 
ties between scale invariant theories to a correspondence 
between theories without scale invariance. Specifically, 
we work with the extension of AdS/CFT to AdS/QCD, 
which has a warped background that breaks conformality 
and gives anomalous dimensions to the operators. The 
scaling dimensions are expressed by a path integral treat- 
ment of propagators of the dual fields in the bulk. Path 
integrals in curved spaces are subtle and we calibrate 
ours using what is known in the AdS/CFT correspon- 
dence. Afterwards we apply the technique to a specific 
warped but asymptotically- AdS background often used 
in AdS/QCD calculations. Second, we calculate the two- 
point function of the gluon field strength operator and 
compare with results computed in quenched lattice QCD. 
In order to carry out our calculation, we must use the 
formula for the flow of the field strength's anomalous di- 
mension into the infrared. In turn, the good agreement 
between our calculation of the propagator and what is 
calculated on the lattice lends strong support to the for- 
mula for the scaling dimension. 

The AdS/CFT correspondence between N = 4 super- 
symmetric Yang-Mills theory in four dimensions and su- 
pergravity on the background AdSs x S 5 , has opened up a 
new venue for probing strongly coupled behavior of quan- 
tum field theories [1]. Since its appearance, there has 
been hope that closely related techniques can be used to 
study strongly coupled physics in theories without super- 
symmetry or conformal symmetry. In particular, there 
has been a great deal of effort to find a gravitational the- 
ory dual to QCD in order to study its low energy, non- 
perturbative regime. Some attempts to achieve this goal 
begin, as in the AdS/CFT correspondence itself, with 
a string theory in a background configuration of branes 
[2, 3]. Others instead posit a gravitational background 
presumed to encapsulate the scaling behavior of QCD. 
The validity of this approach is then tested by using the 
model to compute various quantities in the gravitational 
theory which are compared to experiment results or other 



calculations, often performed using lattice QCD. Agree- 
ment between the two is taken as support for the original 
assumption. This approach to studying low energy QCD 
is known as AdS/QCD [4-7]. In this paper, it is the 
framework we shall adopt. One shortcoming of this ap- 
proach is that it does not address the dynamical origin 
of the backgrounds considered. 

One of the earliest results of the AdS/CFT correspon- 
dence was a computation of the anomalous scaling dimen- 
sions of boundary operators [8, 9]. This was achieved by 
computing the Green functions of bulk fields dual to the 
operators. Let us generalize slightly to a <i-dimensional 
boundary theory and focus on R-symmetry singlet oper- 
ators, meaning we ignore the S 5 coordinates in the bulk 
spacetime. Now the dual gravitational theory effectively 
lives on the background AdS^+i. We use coordinates in 
which the metric takes the form 

ds 2 = g mn dx m dx n = — (dz 2 + 5^ dx»dx v ) . (1) 

Note that we have Euclideanized the metric. A (d — p)- 
form operator on the boundary is dual to a p-form field in 
the bulk. Absent any background fields to provide spe- 
cial reference directions, the index structure of the Green 
function must be trivial. We will write Gab{x,v) — 
$ABG p (x,y), where A and B are generic indices appro- 
priate for a p-form and 

G p (x,y) = (x\(a (p) +m 2 )- 1 \y), (2) 

where D( p ) is the Hodge Laplacian on a p-form field 
[10, 11]. As we will see later, the presence of Wilson lines 
provides a background, and it is only when we neglect its 
background that we keep our trivial index structure. In 
Ref. [8, 9], computation of G p (x 7 y) hinged on the ability 
to solve the relevant differential equations directly, which 
is not possible in the confining background we will con- 
sider later. Instead, we will use a semiclassical evaluation 
of the path integral [12-14]. The Green function is given 

by 
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Here z m (A) is a parametrization of the path followed 
by the particle, e is an cinbcin field associated with the 
parametrization, d + 1 is the dimensionality of the space- 
time and JZ = —d(d+l)/R 2 is the Ricci scalar curvature. 
As emphasized in Ref. [15], there is an ordering ambiguity 
in going from the classical action to the quantum action 
in a curved background. In a spacctime covariant for- 
mulation, the ambiguity allows a curvature term in the 
quantum action with a coefficient that can only be be de- 
termined by imposing some other condition. By dimen- 
sion counting and analyticity, the curvature term must 
be the Ricci scalar. Its coefficient Kd, p is a dimensionless 
constant which can depend only on the spacetime dimen- 
sion and degree p of the form field. We will determine the 
value of Kd,p by reproducing the classical propagator in 
AdSd+i- In the warped but asymptotically AdS^+i back- 
ground, the same value is found by matching results for 
the propagator between two very close boundary points 



because the worldline stays in the near boundary region. 
This value then holds everywhere. 

We make the small fluctuation approximation by ex- 
panding Eq. (3) around its classical solution and treating 
the fluctuations only to quadratic order [14] , which yields 

G{x, y) = D VM (x, y) 1/2 exp[-Sci(x, y)] . (4) 

S c \{x, y) is the classical action of the particle propagating 
from x to y and Dym(x, y) is the Van Vleck-Morette de- 
terminant (VMD), whose square root is the fluctuation 
determinant [14, 16]. The path integral over z depends 
only on s = J Q e dX, so that the path integral over the 
einbein reduces to a simple integral over s. Let z c \ denote 
the classical solution of the equations of motion connect- 
ing x to y and let <r(x,y) be the proper distance of this 
path. Then the classical action term in Eq. (4) becomes 
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for TZmms the Riemann curvature tensor and M mn — 

Thrums 

KiKv The ac tion for the fluctuation field £ comes 
from the second variation of the path length around the 
classical solution. The equations of motion of ( arc 
known as the geodesic deviation equation and the Jacobi 
equation. The particle's trajectory has been decomposed 
as 

z m (A) = z c 7(A) + C m (A), (7) 

with the constraint g m n{z c i)z™( n — so that fluctua- 
tions which are always orthogonal classical solution. This 
means the fluctuation field has only d degrees of freedom. 



AdS<j+i has Riemann curvature tensor 

^mms ~j^2 {dmndrs 9ms9nr^ • (8) 

Taking A to be an affine parameter along the geodesic 
z™(\), the tangent vector i™(A) has constant length, 

g mn (z cl )z%(X)zZ(X) = a(x,y) 2 . (9) 

Using this fact and the orthogonality relationship 
g m n( m Zci = 0, we find 

M mn = -~g mn cr{x, y) 2 /R 2 . (10) 
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Eq. (6) can be expressed using functional determinants 
and computed exactly, 



-DvmO,?/) 1 



1/2 = f Pet det '(-g mn d\ - M mn ) 
Dot det'(-g mn dl) 



<r(x,y)/R \ d/2 
h(a(x,y)/R)J 



-1/2 



(11) 



Here det ' denotes the determinant over vector indices m 
and n only in the subspace orthogonal to the classical 
solution. This means det ' excludes a zero mode from the 
computation of the determinant of M mn . Det denotes 
the functional determinant. There are numerous ways to 
compute the functional determinant but we will evaluate 
it using the Gel'fand-Yaglom Theorem [17, 18], because 
it is used in a less trivial case later. First expand 



Det(-g mn dl 



g mn v(x,y) 2 /R 2 ) 



Det(-g mn dl) 
_ ^Det(-<9; 
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^yf/R 2 ) 



Det(-c^) 

using the diagonality of the metric and cancelled a com- 
mon normalization between the numerator and denomi- 
nator. Even though the theory lives in d + 1 dimensions 
the exponent is only d, reflecting the fact that only fluc- 
tuations tangent to the classical solution are considered. 
The Gel'fand-Yaglom Theorem states that the normal- 
ized determinant of an operator —d 2 — V(X) acting on 
space of functions supported on the interval < A < 1 
with vanishing boundary values is 



Dct(-ff-y(A)) 
Det(-^) 



= «(!)" 



where 



(d 2 x + V(\))u(\)=0, 



(13) 



(14) 



with boundary conditions u(0) = and u(0) = 1. 
Eq. (11) follows immediately for the trivial case of a con- 
stant V = -cr(x,y) 2 /R 2 . 

Using Eqs. (5) and (11) in Eq. (4) shows that the semi- 
classical approximation to the Green function takes the 
form 



G(x,y) = AfK {d _ 1)/2 ia(x,y)^Jm 2 - K d , p TZ 



a(x,y)/R \ d/2 ( ^Jm 2 - K d , P Tl 



(15) 



(d-l)/2 



sinh(a(x,y)/R)J \ a(x,y) 

Af is a normalization factor and Kt d-i)/2 is a modified 
Bessel function. To compute the boundary-to-boundary 
propagator, put x = (x^ 1 , e) and y = (y^, e), for e positive 
but infinitesimal. The geodesic distance is 



<t(x, y) = R In 



\x» - y^\ 2 



+ Q(s) 



(16) 



which diverges as e tends to 0. This pushes the Bessel 
and sinh functions into their asymptotic regimes, with 
sinha; w e x /2 and K( d _iy 2 (x) ~ fix e~ x for any 
d. The asymptotic behavior of the Bessel function, 
the normalization factor in exp(— 5 c i) and the fluctua- 
tion determinant combine to eliminate dependence on 
a(x, y) everywhere except in the exponent. We regularize 
<r{x,y) at both endpoints with some scale /x by setting 
<r(x,y) — CT rog (x,y) — 2i?ln(/ie) and absorbing the di- 
vergent term into the normalization, Af — > The 
resulting Green function is 



Greg(x, y) 



-Af' cxp(-Aa Icg (x,y)/R) (17) 
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The first term on the right comes from Dyyi(x, y) 1 ^ 2 and 
the second term comes from exp[— S c \(x, y)]. By taking 



K d, p 



did ■ 1 : \ 2 l ' 



(19) 



we recover the scaling behavior familiar from AdS/CFT 
[8, 9]. We can make the decomposition 



d(d+l) 



-n^pTl = ^ p(d-p)- 



d - 1 ^ R 
^-Jdjd^l)- (20) 



The first term on the right is carried over from the clas- 
sical field Lagrangian for p-form fields, and comes from 
the Hodge Laplacian [10, 11]. The second term is the 
conformal coupling to the curvature. The third term is 
the path integral measure Jacobian \R~ 2 [19]. If the 
background were instead a sphere of radius R, this term 
would be — \Rr 2 [20]. In general, it adjusts for the dif- 
ference in normalization of the path measure between flat 
and curved spacetimes. 

In Rcf. [21] the Green function is also computed in the 
scmiclassical approximation. The author notes how the 
two-point function of a field in AdS behaves as though 
A/R appears in place of the mass m in the field's clas- 
sical action. Here, we have shown how a path integral 
computation of the two-point function can generate the 
needed terms to convert m into the correct full expres- 
sion for A/R. The path integral technique gives a robust 
understanding of how A emerges from kinematics in a 
more general curved background, which will allow us to 
understand how its value is modified in a warped back- 
ground. 

We now introduce some machinery necessary for the 
computation of the gluon field strength two-point func- 
tion. Yang-Mills field strength operators are not gauge 
invariant, and the field strength two-point function re- 
quires a Wilson loop to form a gauge invariant quantity. 
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We define a loop current on a closed contour in spacetime 
C, parametrized by y(X), by 

j£(x) = J c r(\)SW(x-y(\))d\. (21) 

The current enters the functional integral via the factor 
exp(— 5; n t[A, j]) = exp(— J A ■ j dx), the standard mini- 
mal coupling of a conserved current to a gauge field. In 
non-Abelian theories, the exponential is understood to 
be path-ordered. We use a shorthand notation for the 
expectation value of an operator O in the presence of the 
loop current on C, 

((0))c = N- 1 J [dA] e- s ™ trV{0 exp(- § c A(x)-dx)} . 

(22) 

The trace is over gauge degrees of freedom and V denotes 
path ordering. The normalization is chosen so that the 
value of the partition function Z[jc) = ((l))c — !> 1 as C 
tends to a point. Note that the Wilson loop defined on 
the contour C is precisely Z\jc]- The two-point function 
of field strength operators is 

G$%(x,y\C) = ((F,Ax)F pa (y)))c/((l))c ■ (23) 

Its value depends on C and gauge invariance requires that 
the points x and y must lie on C. Here and throughout 
the rest of the paper, and v refer to boundary direc- 
tions, i.e. all directions except the radial direction z. 

First we review how to compute ((l))c- The non- 
Abelian Stokes' Theorem [221 states 



[dA]: 



-SVm 



Pexp ( - § c A{z) ■ dz) 



[dA] e~ SrM Texp (- § F ■ d 2 a) 



(24) 



Here C is the boundary of surface S with coordinates 
£ and differential area 2-form d 2 a /il/ (£), and F ■ d 2 a = 
\F flv d 2 a tl ' / {^). This latter expression is discretized and 
computed in lattice QCD. On the lattice, path integrals 
over different configurations of the gauge field correspond 
to sums over states with different plaquettes turned on. 
With no dynamical sources of charge around, the plaque- 
ttes which are turned on must span a surface S ending 
on contour C. In the strong coupling limit, the minimal 
area surface dominates [23] . 

The computation of ((l))c can also be performed in 
AdS/CFT [24] and using AdS/QCD [6]. The strongly 
coupled gauge dynamics in the Wilson loop is dual to a 
Nambu-Goto string with endpoints on the contour of the 
Wilson loop. The string worldsheet A m (£) has boundary 
C on the conformal boundary of the (d + l)-dimensional 
spacetime, and 

({l)) c =exp(-S NG [X}). (25) 
Recall the Nambu-Goto action is given by 

S NG [A] = ^7 J d 2 £ Vdet g mn (X)d a X m d b X n . (26) 






FIG. 1. Two possible contours C. The contour on the right is 
the one used in our computations. 



We work in the asymptotically Euclidean AdS^+i "metric 
wall" background of Ref. [6] , with metric 



ds 1 



Qmndx dx 6 



R 2 



{dz 2 + <J dx"dx v ) . 



(27) 

In Ref. [25], this background is compared with alterna- 
tive metrics in its ability to compute rectangular Wilson 
loops, and is shown to provide the best agreement with 
lattice computations of the same. We determine our pa- 
rameters by reproducing lattice results for the Cornell po- 
tential between a heavy quark-antiquark pair modeled as 
a rectangular line current. Doing so yields A ss 330 McV 
and the dimensionless string tension r = R 2 j2-Kol » 
0.1836. The numerator of Eq. (23) is the expectation 
value of two field strength operators inserted along the 
Wilson loop along the contour C in the background gen- 
erated by the rest of C. The insertion of F pl/ at a point 
along C can be achieved by taking the area derivative 
of the contour at this point [26, 27]. In essence, this 
amounts to varying the contour C not with the typical 
Dirac 8 function spike but with a "keyboard variation." 
On the lattice, this variation is made by adding one extra 
plaquette in the fi-v plane at the appropriate point along 
the contour, modifying the contour C to C + SC. We can 
expand the value of the Wilson loop in powers of this 
small change, which gives at first order 

((l))c+8c = Z[ ]c + Sr(x)] a ((1 - a^F^{x)))c ■ (28) 

Here is the oriented area of the extra plaquette added 
at x. Two derivatives insert two field strength operators 



{{F, u {x)F p „{y))) c = 



5 2 Z[j] 



6j^(x)6jp°(y) 



■ (29) 

j=jc 

conserved currents on the 

of gauge fields in the 



In gauge/gravity dualities 
boundary are source terms 
bulk. The keyboard-type variational second derivative in 
Eq. (29) is the two-point function of the field strength in 
the gauge field background with source supported along 
the contour C. It is shown in Ref. [28] for the AdS/CFT 
correspondence that when C is a straight infinite line 

4 



((F^(x)F p(7 (y)))c^l/\x-y\ 4 



(30) 



as expected by the conformality of M = 4 supersym- 
metric Yang-Mills. For theories without conformal in- 
variance, power counting alone does not determine the 
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two-point function. We extend the treatment of the 
field strength two-point function presented in Ref. [28] 
to nonconformal backgrounds. The two-point function 
of a (d — p)-form operator on the boundary corresponds 
to the propagation of a particle transforming as a p-form, 
which we take to have mass to, in the dual gravitational 
theory. In general, form fields have gauge symmetries so 
we will assume for generality that the particle moves in 



a background gauge field A. Its Green function is [12] 

G AB (x,y\A) (31) 
= ( X \{(D{A} 2 + m 2)l + F{A]-n A 1 B \y), 

where £ m " is the generator of spacetime rotations in the 
TO-n plane in the p-form representation, which we take 
to have generic indices A, B. D[A] and F[A] are the 
covariant derivative and field strength computed using 
the background field A only. We can write the Green 
function in terms of a path integral [12], 



z(l)=x 



G AB (x,y\A) = J [<te(A)de(A)]exp(-^ 



e 1 z m z n g mn (z) + e(m 2 - K djt K(z)) 



dX 



(32) 



z(0)=y 



xPexp - / A m {z)z m d\ 



"Pexp ( - ^ F ( z ) • S e rfA 



AB 



The field strength two-point function needs the Wil- 
son loop to be gauge-invariant. In the bulk theory, the 
Wilson loop is dual to a Nambu-Goto string. There- 
fore, gauge invariance dictates the worldline of the field 
strength propagator must lie within the Nambu-Goto 
worldsheet. 

The parallel transporter Pexp(— J A ■ dz) in the field 
strength propagator in Eq. (3) transforms in the adjoint 
representation of the gauge group. It is realized in lattice 
simulation as coincident fundamental and antifundamen- 
tal parallel transporters [29, 30]. The antifundamental 
parallel transporter completes half the line current jc, 
and the fundamental parallel transporter completes the 
other half, as shown in Fig. 1. Thus the worldline of the 
field strength propagator splits the worldsheet into two 
subregions. 

In the computations of Ref. [31], the field strength 
correlator is computed as a weighted sum over the re- 
sults obtained with different Wilson loop contours. For 
this reason, it is not entirely clear exactly which con- 
tour C we should use. Two possibilities are shown in 
Fig. 1. However, there should not be strong dependence 
on the exact contour C in the AdS/QCD calculations, as 
long as the worldsheet it bounds is entirely in the near- 
boundary region of the bulk for small r and can extend 
up to za = (2A) _1 at large r. We use the contour shown 
on the right side of Fig. 1. We take x and y to be sepa- 
rated only in the time direction t by a boundary distance 
of \x — y\ = r. The profile is translationally invariant in 
some boundary direction perpendicular to t. The cross 
section in the t direction is given by z c \{t) which satisfies 
[6] 

„4A 2 Zcl (t) 2 p 4A 2 z m 

; = , (33) 

z cl (;)Vi + <iW 2 z m 



where z m is the maximum value the string's profile as- 
sumes, and is determined implicitly by 




The shape of this profile is shown for several values of r 
in Fig. 2. Notice that for very small r, a slight change 
in r results in a new profile which is just a rescaled ver- 
sion of the old profile. This is due to asymptotic confor- 
mality as z approaches 0. At large r, the string sits at 
z = z. m w (2A) for the majority of its extent in the 
boundary directions. Increases in r serve only to extend 
this portion. Such behavior is common to all choices of 
contour C, with differences lying only in minor details of 
the exact profile. 

There are two ways in Eq. (32) where the background 
field with source supported along the contour C enters 

the evaluation of Greg (x,y\C). First, the background 
field enters through the spin-dependent term 

[VeM-jF--Zed\)] AB = 5 AB + --- . (35) 

We compute only the term proportional to S AB . For a 
2-form, the generic index A = [nv\, so we are keeping 
only the term proportional to 

<V"HH = dupSvo- ~ <W<W 7 (36) 

which is called D±(x,y) in Ref. [31]. 

The other appearance of the background field is 
through 

Vex.p(- $ A • dx) =Vcxp(-jF-d 2 a). (37) 

In the bulk theory, this becomes part of the Nambu-Goto 
string action. To account for its influence, we find the ge- 
ometry of the string worldsheet and particle's worldline 
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which minimize their total action subject to the con- 
straint that the worldline stay within the worldsheet. 
We make the approximation that the balance they strike 
is exactly the string's profile. It is worth repeating for 
emphasis that the point particle's trajectory follows the 
Nambu-Goto string's profile and not the minimum dis- 
tance geodesic connecting the endpoints. To see this, 
consider the quantities 



{{F^(x)F piJ {y))) c = exp(-S NG [A(C)] - S ci (x,y)) (38) 
«l» c = exp(-S NG [X(C)]) (39) 



Although the Nambu-Goto action cancels out of 

G$$(x,y\C) = ((F^(x)F pa (y)))c/((l))c in the term 
oc 5[^][ pcr ], it still affects the geometry of the minimal 
action state. It determines it completely in our approx- 
imation. There is interaction energy which will cause it 
to deviate from the case without field strength operator 
insertions along the boundary, but the deviation from 
our approximation is very small. So we must only com- 
pute the action of the field strength propagating along 
the string's profile. 

When we drop the spin-dependent coupling term in 
Eq. (32), the propagator once again takes the form in 
Eq. (4). But now the curvature terms vary with position 



cxp 



and 



r 

-S c i(x,y)] = 

Jo 



ds 



(27r S )( d +l)/ 2 ^ ^ 2S 



a{x,y) 2 s f 1 



-\f [m 2 - K d . p 1l(z cl {X))}dXj 



(40) 



(d-l)/4 



= 2 ( Jot™ 2 _ K d,pR-{z cl (\))]d\l a(x,y) 2 j K (d-i)/2 f <?(x, y)\J [m 2 - K d , p TZ{z cl (X))]dX 



C(i)=o 



DvMix^y) 1 ' 2 ^ I [dC(A)]exp 



C(o)=o 
C(i)=o 
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2s Jo 



9mn(Zcl)C C 7^"mpnq{Zcl)C ^cl^ ^cl 



dX 



(41) 



= J [d((X)}exp^j\ m {g mn d 2 +M mn }C n dXy 



C(o)=o 



The field strength operator is a 2-form, so for d = 4 
there is a simplification because — 0. Moreover, 
its mass vanishes. This can be determined in the near 
boundary region, where A = 2 and the AdS / CFT scaling 
dimension formula applies. However, we must first take 
the large cr(x,y) limit of Eq. (40), keeping m finite as 
a regulator, and afterwards we can absorb m into the 
normalization. At the end of these manipulations, we 
find 

exp[-S cl (x, y)} -+ Afa(x, y)- d ' 2 . (42) 

DvMiXjy) 1 ^ 2 can be expressed using the determinant 
of the operator g mn d\ + M mn , 



Detdet'(- 9mra <9|-M„ 
Dct det'(-fif m „^) 



-1/2 



(43) 



In the background metric of Eq. (1), the curvature is 
constant and isotropic. This means M mn (z c i(A)) is a con- 
stant as in Eq. (11). This feature also holds to a very 
good approximation in the near-boundary region of the 
warped background in Eq. (27). Furthermore, for large 
r = \x — y\ the string's profile will sit for the majority of 
the time at za = (2A) _1 , and for these points the value 



of M mn (z ci (X)) varies extremely slowly. The variation of 
M mn with A in the intermediate region is also slow. The 
curvature tensor in Eq. (51) is bounded above in abso- 
lute value by the unwarped curvature tensor in Eq. (8). 
The curvature is taken to be small compared to the scale 
set by the string tension. So it is justified to expand in 
(d 2 )~ x M mn , and we use the expansion of the determi- 
nant of nearly unit matrix, det (1 +e) = 1 +tre + 0(e 2 ), 
to write [32] 



det'{-g mn dl-M mn ) (44) 
= det'(-<w^) x [l + (dl)-hiM] + 0(M 2 ) 

= (-d\) d - i-dl)^ 1 tr M + 0(M 2 ) 
= (-df - d-HvM) d + 0(M 2 , d 2 M) 



We then drop the higher order terms. For an unwarped 
AdS metric, the above result is exact and the dropped 
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FIG. 2. String profiles with various values of r 



terms in fact vanish. Next, use the relationship 



D VM {x iy ) 1/2 - 



Bct[(-d 2 x -d-HrM mn ) d } \ - 1/2 
~ Det[(-d^] 



(45) 



Det(-^-rf- 1 trM TOra ) 
Bet(-dl) 



-d/2 



This manipulation is valid here and in Eq. (12) because 
the functional determinant is a power of a single differ- 
ential operator and not of the product of different differ- 
ential operators [33]. The functional determinant is now 



readily computed using the Gel'fand-Yaglom theorem, 



Det(-02 -d-HrM mn ) 



Det(-dl) 
where u(X) satisfies 

(dl + d-hrM mn )u(\) =0 



(46) 



(47) 



with boundary conditions u(0) = and ii(0) = 1. Since 
we are taking M mn to be small and slowly varying, we 
use a WKB approximation, resulting in 



£,1/2 



sinh(/o A /-d- 1 trM mn (A)dA Q ^ ~ d/2 
^-d-HrM mn (0) 



(48) 



To express V—d^tr M mn , we first introduce a frame field 
e m a = {R/z)5 m a or equivalently e m a — (z/R)S m a , so 
that g mn = e m a e n b S ab . We write 



z^(X) = a(x,y)e m a t a (X) 



(49) 



where i a (X) is the tangent along the geodesic normalized 
in the local rest frame, 5 a b i a (X)i b (X) = 1. For notational 
convenience, we introduce 9(X) which satisfies 

cos 2 <9(A) = t x (A) 2 and sin 2 (9(A) = t z {\) 2 . (50) 



The Riemann tensor of the metric in Eq. (27) has components 



K 



n 



1 

1 



(1 - 4A z ) (g w g va - g^gvp) 



(l + 4A 2 z 2 )g^g, 



(51) 



and all others vanish except those related by symmetries of the Riemann tensor to the two components above. We 
find that 



v /-d- 1 trM mn (z cl (A)) = ^ e- 2A " z2 J (1 - d-^cos 2 9(1 - 4A 2 z 2 ) 2 + sin 2 0(1 + 4A 2 z 2 )] + dr^l + 4A 2 z 2 ) , (52) 

R * 

where z — z^(X) is the radial coordinate at a point along the classical path. 

I 

As x and y approach the boundary, sinh assumes its Using an affine parameter in the expressions above allows 
asymptotic form and we arrive at a simplification like the factorization of the distance cr(x, y) out of the integral 
in Eq. (17), so that the Green function takes the form 

G^ FF \x, y\C) = AAcxp[-A • a(x, y)/R] , (53) 

where A is defined by 



in Eq. (55) . For a generic parametrization of the classical 
solution with parameter £, the quantity 

a-A/R= r ^-d-HrM(z cl (0)\\dz cl (0\\ (56) 



and 



A(z cl (A)) = R^y/-d- 1 tiM mn {z d {X)) (54) 



u ■ A/ R = f A(z cl (A)) dX . (55) 



where ||dz c i(OII = vt^dPiK). z ci (£ x ) = x and 
z c\(^y) — V- We will perform calculations using the t- 
coordinate parametrization of the classical solution given 
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in Eq. (33), for which 

r/2 

«7-A/i2= | «)^^ + <i(t) S *- (57) 

-r/2 

The divergence of the path length as the endpoints ap- 
proach the boundary must be regularized like in Eq. (17), 
so that 

G£p{x,y\C) = W exp[-A • a ieg (x,y)/R] (58) 

Notice that A • <7 rcg is not just a function of the z- 
coordinate of the classical solution but also the tangent 
to the classical solution. A is the scaling dimension of the 
gluon field strength and the z-coordinate is dual to the 
physical extent of field configurations on the boundary. 
This means the anomalous scaling dimension is not only 
scale-dependent but is sensitive to rate of change of the 
physical scale as the particle propagates. 

When using a generalized proper time regularization of 
multidimensional path integrals, care must be taken to 
handle the multiplicative anomaly, which is the difference 
between results of the two orders in which one might take 
the functional and finite determinants in Eq. (43) [34, 35]. 
Some authors [32, 36] claim either explicitly or implicitly 
that it is correct to take the functional determinant of 
g m nd\ + M mn before taking the determinant over finite 
indices, i.e. computing 



/ det ' Det{-g mn d\ - M mn ) 



det'Det(~g mn dl) 



-1/2 



(59) 



instead of Eq. (43). To evaluate Eq. (59), one can use a 
multidimensional Gel'fand-Yaglom Theorem [32], 



det ' Det(-g mn d| - M mn ) 
det ' Det(-# mjl c^) 

where A m n (X) solves 



det^l)- 1 



(g mn d 2 x + M mn )A n p (\) = 



(60) 



(61) 



with boundary conditions A m n (0) = and A n (0) = 
S m n - We will not pursue the details, but using a WKB 
approximation of this differential equation results in a 
propagator of the form in Eq. (58) except with 



A(z cl ) = i J Rtr v /-M m „(z cl ) 



(62) 



instead of Eq. (54) . Ref. [35] shows that the order used in 
the present paper is the correct order, and gives results in 
agreement with other non-functional methods. Further- 
more, we will see that the use of Eq. (62) does not give 
a propagator in good agreement with the lattice data. 

We can now compare the accuracy of the AdS/QCD 
calculations with lattice results. We arrive at an explicit 
expression for Greg^ (x, y\C) when we combine Eqs. (52), 
(54) and (58). For C we are using the contour shown on 



Log 10 [G FF (r)] 
11 

10 

9 
8 
7 




0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



r/fm 



FIG. 3. Boundary-to-boundary field strength correlator with 
the correct A (solid) from Eq. (54) contrasted with the in- 
correct A (dashed) of Eq. (62), and A = 2 held constant 
(dotted). Lattice data was taken from Ref. [31] and error 
bars from Ref. [37]. 



the right in Fig. 1, which produces z c \ given in Eq. (33). 
The advantage of this choice is that we can write explic- 
itly 



cos 2 9(t) 



l + <i(*) 2 

{z c i(t)/z m ) 4 exp[8A 2 ;4 



(63) 



8A 2 z cl (i) 2 



Fig. 3 shows the AdS/QCD computation of the propaga- 
tor with d = 4 plotted against the lattice data [31, 37]. 
The solid curves show the propagator Eq. (58). For con- 
trast, if we use the incorrect ordering of determinants to 
compute A as in Eq. (62), the dashed curve of Fig. 3 
results. At small r, the string sits entirely in the near 
boundary region, as seen in Fig. 2. There, the asymp- 
totically conformal behavior dictates G(r) cx l/r 2Auv , 
where A uv = A(z — > 0) = 2 for the field strength oper- 
ator. The dotted curve in Fig. 3 shows the propagator 
computed using Eq. (57) but with the scaling dimension 
fixed everywhere at the ultraviolet value, A(z c i(A)) = 2. 

The normalizations of the dashed and dotted curves 
were set to agree with the lattice data in the small r 
limit, where the universal 1/r 4 behavior must hold. For 
the solid curve, the normalization is determined by min- 
imizing the x 2 value of the fit. The end result gives 
X 2 /d.o.f. w 3.07. 

At larger distances, starting around r = 0.4 fm, the 
warp factor has a notable effect on the string's profile. In 
unwarped AdS and in the small z regime of the metric in 
Eq. (27), string profiles look simply like rescaled versions 
of string profiles with the endpoints closer together. Such 
behavior ceases to hold for large enough r. In fact, soon 
the string nearly saturates the bound z m < za = (2A) _1 . 
Further increases in r simply extend the portion of the 
string's profile which sits at the maximum z value. This 
yields the exponential behavior of the propagator at these 
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distances, 

G(a:,2/)ocexp(-|a:-2/|/Agi), (64) 
where the gluonic correlation length is 

A g i = R (A(z m ) v /ffoo(z m ) y 1 = (2V2A)- 1 . (65) 

Substituting A w 330 MeV, we find that 

A g i w 0.21 fm . (66) 

Our result is in good agreement with the lattice results 
of Ref. [31], which indicate A g i w 0.22 fm. Ref. [38] has 
A w 0.11 — 0.13 fm, different roughly by a factor of 2. 

In this paper, we have addressed a few closely related 
questions. We set out to compute the two-point func- 
tion of the gluon field strength operator in the absence 
of dynamical quarks or other sources of color charge in 
the fundamental representation. The result is known 
from quenched lattice QCD computation. To perform 
this calculation, we first had to gain an understanding of 
the origin of the scaling dimension formula known from 



AdS/CFT. From there, we learn how to modify it ap- 
propriately for a theory which flows with the physical 
scale. The results presented here should open up a num- 
ber of lines of further inquiry, including how to extend 
to more complicated operators, higher order correlation 
functions, or making systematic improvements to the ap- 
proximations contained in this paper. 
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